Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN MfeK style file vl.4) 



Studies of CMB structure at Dec=+40°. 
II: Analysis and cosmological interpretation 

S. Hancock 1 , CM. Gutierrez 2 , R.D. Davies 3 , A.N. Lasenby 1 , G. Rocha 1 , 
R. Rebolo 2 , R.A. Watson 2 ' 3 and M. Tegmark 4 

1 Milliard Radio Astronomy Observatory, Cavendish Laboratory, Madingley Road, Cambridge CB3 OHE, UK 

2 Instituto de Astrofisica de Canarias, 38200 La Laguna, Tenerife, Spain 

3 University of Manchester, Nuffield Radio Astronomy Laboratories, Jodrell Bank, Macclesfield, Cheshire, SK11 9DL, UK 
4 Max- Planck- Institut fur Physik, Fohringer Ring, D-80805 Miinchen, Germany 



ABSTRACT 

We present a detailed analysis of the cosmic microwave background structure in 
the Tenerife Dec=+40° data. The effect of local atmospheric contributions on the de- 
rived fluctuation amplitude is considered, resulting in an improved separation of the 
intrinsic CMB signal from noise. Our analysis demonstrates the existence of common 
structure in independent data scans at 15 and 33 GHz. For the case of fluctuations 
described by a Gaussian auto-correlation function, a likelihood analysis of our com- 
bined results at 15 and 33 GHz implies an intrinsic rms fluctuation level of 481^ /iK 
on a coherence scale of 4°; the equivalent analysis for a Harrison-Zel'dovich model 
gives a power spectrum normalisation of Qrms-ps — 22^g° /iK. The fluctuation am- 
plitude is seen to be consistent at the 68 % confidence level with that reported for the 
COBE two-year data for primordial fluctuations described by a power law model with 
a spectral index in the range 1.0 < n < 1.6. This limit favours the large scale CMB 
anisotropy being dominated by scalar fluctuations rather than tensor modes from a 
gravitational wave background. The large scale Tenerife and COBE results are consid- 
ered in conjunction with observational results from medium scale experiments in order 
to place improved limits on the fluctuation spectral index; we find n = 1.10 ± 0.10 
assuming standard CDM with H = 50kms _1 Mpc \ 

Key words: Cosmology - Large Scale Structure of the Universe - Cosmic Microwave 
Background 



1 INTRODUCTION 

Observations of fluctuations in the Cosmic Microwave Back- 
ground (CMB) have been widely recognized to be of funda- 
mental significance to cosmology, offering a unique insight 
into the physical conditions in the early Universe. The am- 
plitudes and distribution of such fluctuations provide criti- 
cal tests of the origin of the initial perturbations from which 
the structures seen today have formed. On scales > few de- 
grees, CMB observations probe scales of 1000's of Mpc, in- 
accessible to conventional astronomy. At these large angles, 
the structures form part of an intrinsic spectrum of fluctu- 
ations generated through topological defects or inflation. In 
this linear growth regime, observations of the scalar CMB 
fluctuations provide a clean measure of the normalisation 
of the intrinsic fluctuation power spectrum. This normal- 
isation has been established by a number of independent 
CMB observations (Smoot et al. 1992, Ganga et al. 1993, 



Hancock et al. 1994). In many theories, tensor CMB fluc- 
tuations from a background of gravitational waves are also 
expected to be significant on these large scales, and mea- 
suring the slope of the power spectrum offers the potential 
to constrain this contribution to the CMB anisotropy (Han- 
cock et al. 1994, Steinhardt 1993, Crittenden et al. 1993). A 
comparison of the large-scale anisotropy results with those 
on medium scales can also be used to separate the scalar 
and tensor components under the assumption of a specific 
cosmological model. 

The Tenerife CMB experiments were initiated in 1984, 
with the installation of the first 10 GHz switched-beam ra- 
diometer system at the Teide Observatory on Tenerife Is- 
land. A subsequent programme of development has led to 
the present trio of independent instruments working at 10, 
15 and 33 GHz. The ultimate objective is to obtain three 
fully sampled sky maps covering some ~ 5000 square de- 
grees of the sky at each frequency and attaining a sensitiv- 
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ity of ~ 50 jjK at 10 GHz, and ~ 20 /iK in the two highest 
frequency channels. Drift scan observations have been con- 
ducted over a number of years covering the sky area between 
Dec=+30° and +45°. The deepest integrations have been 
conducted in the Dec=+40° region and resulted in strong 
evidence for the presence of individual CMB features (Han- 
cock et al. 1994). 

Davies et al. 1995 (hereafter Paper I) described the per- 
formance of the experiments and gave an assessment of the 
atmospheric and foreground contributions to our data at 
Dec=+40°; here we analyse in detail the results and cos- 
mological implications of such observations. Section 2 de- 
scribes the observational strategy and presents the stacked 
scans at each frequency. In Section 3 we use several statisti- 
cal methods to calculate the level of the detected signals and 
their origin. A statistical comparison with the results of the 
COBE DMR two-year data is conducted in Section 4 and 
used to place limits on the spectral index of the primordial 
fluctuations. In Section 5 we use the additional information 
provided by medium-scale anisotropy results to provide im- 
proved limits on n. 



2 THE SCANS AT DEC +40° 
2.1 Observations 

Observations were conducted at the three frequencies 10, 15 
and 33 GHz by drift scanning in right ascension at a fixed 
declination of 40°. The measurements were made indepen- 
dently at each frequency, using separate dual-beam radiome- 
ter systems as described in Paper I. The three instruments 
are physically scaled so as to produce approximately the 
same beam pattern (FWHM~ 5°) on the sky, thus allow- 
ing a direct comparison of structure between frequencies. 
A characteristic triple beam profile (switching angle 8°1) is 
obtained by the combination of fast switching (63 Hz) of 
the horns between two independent receivers and secondary 
switching (0.125 Hz) provided by a wagging mirror. We make 
repeated observations of the sky, binning the data in 1° in- 
tervals in RA and stacking them together in order to reduce 
the noise as compared with individual measurements. As a 
consequence of using two independent channels, the receiver 
noise contribution to the final data scans is reduced by a 
factor y/2 compared with single channel observations. 

The data considered in this paper are the same as 
those presented in the preliminary report by Hancock et al. 
(1994). As before we restrict our analysis to the RA range 
161° — 230° corresponding to Galactic latitude b > 56°. At 
these high latitudes foreground emission from the Galaxy is 
expected to be at a minimum. This sky region has also been 
selected (see Paper I) to be free from discrete radio sources 
above the 1.5 Jy level at 10 GHz; assuming a flat spectrum 
this corresponds to an antenna temperature contribution of 
48, 21 and 5 fiK at 10, 15 and 33 GHz respectively in the 
Dec=+40° data. The actual contamination is reduced be- 
low these levels by using the Kiihr catalogue (Kiihr et al. 
1981) supplemented by the VLA calibrators list to subtract 
the discrete source contribution from each data scan. This 
amounts to making corrections to the rms signal level of 
13, 8, and 3 ^iK at 10, 15 and 33 GHz respectively The 
contribution of unresolved radio sources is expected to be 




220 



15 GHz 



200 180 

RA (degrees) 



160 



a. 



tn 
CD 
o 
C 
<p 
u 

CD 



-a 
a 

o 
o 

CD 

m 




220 



33 GHz 



200 180 

RA (degrees) 



160 




220 



200 180 

RA (degrees) 



160 



(15+33) 




220 200 180 

RA (degrees) 



160 



Figure 1. The stacked scans at Dec=+40°. The second differ- 
ences in temperature arc shown binned at 4° intervals in RA. The 
error bars are at 68 % confidence and take into account any atmo- 
spheric correlations between receiver channels. The scans labeled 
15+33 represent the weighted addition of the data at 15 and 33 
GHz. 



AT/T < 5 x 10~ 6 at 15 GHz and significantly smaller at 33 
GHz (Aizu et al. 1987, Franceschini et al. 1989). 

Figure 1 presents our final stacked scans at 10, 15, 33 
GHz and the weighted addition (named 15+33) of the data 
at the two higher frequencies. We have calculated the stan- 
dard error-bars of each point considering the presence of cor- 
related atmospheric noise between the two channels of each 
receiver system as described in Paper I. This represents a 
mean enhancement of 3 % at 10 GHz, 4 % at 15 GHz, and 
16 % at 33 GHz with respect to the estimation of the noise 
computed in Hancock et al. (1994). The revised sensitivities 
per beam-sized area are 61, 32, 25 and 20 fiK at 10, 15, 
33 and 15+33 respectively. A visual inspection reveals the 
presence of common features in the scans at the two higher 
frequencies; this is supported by a cross-correlation analysis 
in Section 3.2. 
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2.2 Reliability of the detected signals 

Determination of the amplitude of the CMB component of 
the structure requires one to consider the contributions of 
random noise and foreground signals to the observed data 
scans. The former has its origin in the thermal variations in 
the receivers and in the fluctuating component of the atmo- 
sphere, whilst the latter consists primarily of free-free and 
synchrotron emission in the Galaxy, plus emission from the 
Sun and Moon. Of these effects, only the Galactic emission 
remains constant from day to day at a given frequency. In 
Paper I we have estimated a maximum Galactic contribu- 
tion of AT rms — 4/iK in the results at 33 GHz. An improved 
separation between the Galactic and the cosmological signal 
at each frequency will be presented in a forthcoming paper. 
We have removed the effects of the Sun and Moon to better 
than 50 dB in each daily scan, and stacking n days of ob- 
servations with the Sun and Moon in different positions will 
reduce any residual contribution to an insignificant level. 

The thermal noise contribution from the receivers con- 
forms to a random Gaussian distribution and integrates 
down accordingly. In each instrument, the channel 1 receiver 
is fully independent of the channel 2 receiver and hence 
there is no correlation between the receiver noise in the data 
recorded for the two channels. Atmospheric emission varies 
randomly over days and thus at each frequency, stacking to- 
gether n days of observations reduces the atmospheric sig- 
nals in a given day by ^Jn. Additionally, the Tenerife observ- 
ing strategy is such that the independent 10, 15 and 33 GHz 
instruments observe different declination strips on any given 
day so that the atmospheric contribution to data taken at 
each frequency is uncorrelated. This offers a key advantage 
over experiments such as Saskatoon (Wollack et al. 1993), 
where the data in all of the frequency channels are highly 
correlated thus introducing an element of uncertainty due 
to the need to account for this. In our individual measure- 
ments the atmospheric contribution to each channel is cor- 
related, since the timescale of the receiver switching is short 
compared with that for typical atmospheric fluctuations. In 
previous papers (Davies et al. 1987, Watson et al. 1992, Han- 
cock et al. 1994 and Gutierrez et al. 1995) this effect was not 
considered. The presence of a correlated atmospheric compo- 
nent affects the estimation of the astronomical signal based 
on the difference of the (A + B)/2 and (A — B)/2 data sets 
(A and B being separate subsets of the data). In particular 
if A and B each correspond to one channel, (A + B)/2 con- 
tains the contribution of the astronomical signal, and the 
atmospheric and instrumental noise, whilst (A — B) /2 only 
contains the instrumental noise. This is the case for the anal- 
ysis of the 33 GHz data presented in Hancock et al. (1994) 
and therefore the signal level obtained from the difference 
in variance between the (A + B)/2 and (A — B)/2 scans 
contains a contribution from the atmosphere as well as the 
astronomy. Note that at 15 GHz the data were subdivided 
according to the observing epoch and consequently the con- 
tribution of the atmosphere to the derived signal was largely 
reduced. 

Here we present a new split of the 33 GHz data into 
two subsets X and Y such that both channels of a given 
scan are included in the same data subset. Considering the 
non-repeatiblity of the atmosphere from day to day com- 
mon atmospheric signals are not expected to occur in both 



subsets. Also, considering the large number of independent 
observations (more than 50 for each channel) we expect that 
the net effect of the atmospheric signals will be just an in- 
crease of the variance in the final stacked scans of each sub- 
set, this increase in variance being approximately the same 
in both. For these reasons, when X and Y are combined to 
form the (X + Y) /2 and the (X — Y) /2 scans we expect 
to have approximately the same atmospheric signals in the 
sum and in the difference. Figure 2 shows the stacked data 
scans for the X and Y subsets and their sum and difference. 
We see a general agreement between the results of splits X 
and Y , and the presence of common features in both. An 
analysis similar to that presented in Hancock et al. (1994) 
gives an astronomical signal (a s — <J^x+y)/2 — cr (x-Y)/2) 
with an amplitude a s = 43 ± 12 /iK. The value of the sig- 
nal quoted in Hancock et al. was a id = 49 ± 10 fiK; the 
difference with our improved estimation is certainly due to 
the subtraction of the atmospheric signal. The difference be- 
tween both estimates (a = (p'o ld — al ) 1//2 ) is ~ 23 /iK which 
is our best assessment of the atmospheric contamination in 
the analysis based on the addition and difference of the 33 
GHz data; this value is in good agreement with estimates 
obtained using other methods. 



3 STATISTICAL ANALYSIS 

The estimates of the astronomical signals made in the pre- 
vious section can be improved on using a detailed statistical 
analysis. Here we use a likelihood analysis and take into ac- 
count the contribution of the correlated atmospheric noise 
by enlarging the error bars on the scans as shown in Figure 
1. These modified scans are also used to study the presence 
of common features at 15 and 33 GHz by the calculation of 
the cross-correlation function. 



3.1 Likelihood analysis 

This analysis takes into account all the relevant parameters 
of the observations: experimental configuration, sampling, 
binning, etc. The combination of the atmospheric and in- 
strumental noise can be modeled by a Gaussian distribu- 
tion uncorrelated from point to point, (see previous section), 
which implies that the noise only contributes to the diago- 
nal terms of the covariance matrix. We also assume that the 
astronomical signal is described by a Gaussian random field 
and therefore our results correspond to a superposition of 
Gaussian fields in which all their statistical properties are 
specified by the covariance matrix, which takes into account 
the full correlation between the data points. 

We have analyzed our results for two hypothetical sky 
models, the first of which corresponds to a signal described 
by a Gaussian auto-correlation function (ACF) with am- 
plitude \/Co and width 6 C . This is not a realistic physi- 
cal scenario but has been used widely in the past (Davies 
et al. 1987, Readhead et al. 1989, Watson et al. 1992) be- 
cause it provides for an easy comparison between the results 
of experiments with different configurations. The intrinsic 
ACF for these models is given by 

C(0) = O,exp(-^). (1) 
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Figure 2. The data scans obtained from the analysis of two sub- 
sets of the 33 GHz data. In a) and b) arc shown the subsets de- 
noted X and Y which arc constructed such that the atmospheric 
contributions to each are independent: see main text for details. 
The differences in the (X + Y)/2 (panel c)) and (X-Y)/2 (panel 
d)) scans should therefore be due only to astronomical signals. 



which is modified accordingly by our triple beam filtering 
(Watson et al. 1992). Our instrument is sensitive over a 
range of coherence angles 1° <9 C <10° , attaining peak sen- 
sitivity for a coherence angle of 4° . We have analyzed the X 
and Y subsets for 15 and 33 GHz, the total stacked scans at 
these two frequencies, and our best scan 15+33. The results 
for 9 C — 4° are presented in the third column of Table 1. 
The amplitude of the intrinsic signal corresponding to the 
maximum likelihood is given, along with the one-sigma con- 
fidence bounds calculated in a Bayesian sense with uniform 
prior. All results look consistent with clear detections at 
the two to three sigma level, and mean values of the signal 
slightly smaller than those presented in Hancock et al. , due 
to our improved estimate of the error bars in the stacked 



scans (see Section 2 and Paper I). Figure 3 presents the 
contours of equal probability for the 15+33 scan. We see a 
well defined point of maximum likelihood at 6 C ~ 4° and 
v/Co ~ 50 fiK. 

The second model considered here is more interesting 
from a cosmological viewpoint. It corresponds to the predic- 
tion of the power law form (P(k) oc k n ) for the spectrum 
of the primordial fluctuations. Considering only the Sachs- 
Wolfe part of the spectrum of the fluctuations the intrinsic 
ACF can be expressed as 

C(8) = Qrms ~ ps2 V(2Z + l)C, (s) fl(co60) (2) 
i 

MS) = (S) r[Z + (n-l)/2]r[(9-n)/2] 
; 2 r[i + (5 - n)/2] r[(3 + n)/2] 

where the sum is extended to the multipoles I < 60 which 
corresponds to the range of angular sensitivity of our ex- 
periments. For I > 20, standard models predict additional 
contributions to the CMB anisotropy, as one moves into the 
low I tail of the CMB Doppler peak. Hence fitting for the 
Sachs- Wolfe term alone (Equation 2) to CMB data on these 
scales can lead to the derived values for n being increased by 
as much as 10% over the true primordial value. This point 
should be borne in mind when comparing the limits on n 
from the Sachs- Wolfe term (Sections 3 and 4) to those from 
a fit to a full CDM type functional form as in Section 5. 

For a given value of the spectral index n, the intrinsic 
ACF is a function only of Q r ms- P s ■ Figure 4 shows the like- 
lihood surface as a function of Q r ms- P s and n for the 15+33 
scan. The peak likelihood forms a ridge displaced from zero 
in Qrms-ps and corresponds to a 3 — 4 sigma detection of 
structure for each value of n considered. The shape of the 
surface implies that all values of n in this range are equally 
likely. This is predominantly a consequence of our observ- 
ing technique which samples only a small angular range of 
the spectrum of fluctuations. Thus whilst our observations 
provide a good measure of the fluctuation power on ~ 4° 
scales, they do not in themselves contain sufficient informa- 
tion about the distribution of power with angular scale to 
allow a useful determination of the spectral slope: for this 
one must compare with experiments on other angular scales 
(see Sections 4 and 5). For the specific case of a Harrison- 
Zel'dovich spectrum (n = 1) the results of the likelihood 
analysis are given in column two of Table 1. The normalisa- 
tion Qrms-ps corresponds to the maximum of the likelihood 
function and the confidence intervals are at 68 %, calculated 
in the standard Bayesian manner using uniform prior. We 
see that in general the results are consistent and agree with a 
global normalization of the quadrupole Q rms _ ps ~ 20 — 25 
1-iK. Our best estimate for Q r ms- P s of 22+ 6 fjK from the 
15+33 scan is reduced over the value of 26 ± 6/iK previously 
reported due to our now having properly accounted for the 
correlated atmospheric noise. 



3.2 Cross-correlation analysis 

We have investigated the presence of common features in 
our 15 + 33 scans at Dec=+40° by evaluating the cross- 
correlation function, 
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Table 1. Results of the likelihood analysis for a Harrison- 
Zcl'dovich spectrum of fluctuations (second column) and for a 
Gaussian ACF (third column). 

v (GHz) Q rms — ps ( M K) (pK) 

ISA 27±{ 6 6 54+J? 

15B 17± 9 17 24±H 

15 2ltl 2 44+26 

33A 22±\* 45+J 2 

33B 28+g 2 57+Jg 

33 24+i 1 49l^ 
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Figure 3. Contour levels of equal likelihood for the 15+33 scan 
in the case of a Gaussian shaped ACF. The contours correspond 
to 10, 20, 30 ... 90 % of the total probability distribution. The 
Tenerife configuration obtains maximum sensitivity for coherence 
angles in the range 2°^9^,6°\ structure is clearly detected at ~ 
50/iK over this angular scale range. 

C(9) = ^= (3) 

where AT and AT' denote the second differences in the 
stacked scans at 15 and 33 GHz respectively. The sum is ex- 
tended over all points subtending a mutual angle 6 in RA, 
with weights w» = 1/af and w'j = 1/af respectively. The 




Figure 4. The two-dimensional normalised likelihood surface as 
a function of the spectral index n and the normalisation Q r ms-ps 
for the Tenerife data. The projected contours are at 68 %, 95 % 
and 99 % confidence. 



data points in Figure 5 represent the mean and 68 % confi- 
dence intervals of C(8). The observed profile is characteristic 
of the triple beam sampling and it is clear that for angles 
8 < 20° the results are consistent with the presence of com- 
mon signals in the two independent data sets. Because of the 
non- independence of the data points in the cross-correlation 
function, the statistical significance is difficult to evaluate, 
requiring the computation of the covariance matrix of the 
correlation function (see e.g. Ganga el al. 1993). However, 
considering only the correlation at zero lag, the observed 
value of C{9) — llOO+S^S' M^ 2 is found to be inconsistent 
with noise at the 95 % level. This value for the amplitude 
of the component of the signal common to the 15 and 33 
GHz data is in agreement with the equivalent estimates of 
the variance of the signals present in each data scan sepa- 
rately, demonstrating that common signals can account for 
all of the structure seen in the separate scans. This is con- 
trary to what one expects for structure of Galactic origin, 
which in the case of synchrotron and free-free emission would 
exhibit a change in amplitude by a factor 9 and 5 respec- 
tively on moving from 15 GHz to the higher frequency of 
33 GHz. This result taken in conjuction with the limits on 
the Galactic contamination of the 33 GHz scan (see above) 
and the successful comparison of the features seen in the 15 
and 33 GHz scans with the higher frequency COBE data 
(Lineweaver et al. , 1995), leads us to assign a cosmological 
origin to the majority of the signal present in the 15 GHz 
data. 

We have explicitly tested the consistency of the com- 
mon signals with fluctuations originating in inflationary 
type models which predict an approximately scale invari- 
ant n ~ 1 spectrum of fluctuations. Using the normalisation 
of Qrms-ps = 22/iK as derived from the likelihood analy- 
sis of the 15+33 scan we simulated 5000 realisations of a 
sky conforming to a Harrison- Zel'dovich spectrum and re- 
peated the Tenerife sampling strategy. The mean value (solid 
line) and one-sigma confidence intervals (dashed lines) of 
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the auto-correlation function obtained from these realisa- 
tions are plotted in Figure 5. There is clearly a good level of 
agreement in the amplitude and shape of the experimental 
results and the theoretical model. Statistically the Tener- 
ife observations have the properties expected for primordial 
fluctuations generated from inflation, and although these 
results do not prove inflationary theory, this does offer a 
plausible explanation for the existence of structure on scales 
greater than the horizon volume at recombination. In this 
case, the clearly defined structures visible in the scans in 
Figure 1 represent the seed perturbations generated in the 
inflationary era at t < lCF 34 s. In the following section we 
investigate the potential to prove or disprove inflation by 
measuring the slope n of this primordial spectrum of fluctu- 
ations using large scale CMB anisotropy measurements. 



4 STATISTICAL COMPARISON WITH COBE 
DMR 

A comparison between the results of different CMB experi- 
ments offers the opportunity to check independent measure- 
ments, to extend the range in frequency and angular scale, 
to constrain cosmological models and, if the sensitivity of 
the experiments is sufficient, to compare features. There are 
several experiments operating on angular scales of a few de- 
grees: MIT (Ganga et al. 1993), COBE (Smoot et al. 1992, 
Bennett et al. 1994), RELIKT (Strukov et al. 1993), ARGO 
(De Bernardis et al. 1994a) and Tenerife. The first compar- 
ison between independent CMB observations was made by 
Ganga et al. (1993) who found a clear correlation between 
the results of the first year of COBE DMR observations and 
those of the MIT experiment. De Bernardis et al. (1994b) 
have also made a statistical comparison between the am- 
plitude of the signal reported for ARGO and that of the 
COBE DMR first year results, from which they constrain 
the spectral index to be 0.5 < n < 1.2 in the absence of any 
gravity wave background. In our preliminary report (Han- 
cock et al. 1994) we compared the amplitude of the signal 
detected on ~ 5° scales in our Dec=+40° data with that 
on ~ 7° scales for the first year of COBE data. We found 
that both results were consistent with inflationary models 
(n>0.9) but with a favoured spectral index of n = 1.6. 
Lineweaver et al. (1995) have presented the first direct com- 
parison of CMB features between the two-year COBE DMR 
data and the Tenerife Dec=+40° observations, confirming 
the agreement in the level of the normalization of both ex- 
periments and providing clear evidence for the presence of 
common hot and cold spots in both data sets. The compari- 
son presented here is different to that conducted in Hancock 
et al. (1994) in that we use a more rigorous comparison 
technique that utilises the likelihood function to incorpo- 
rate fully the effects of cosmic and sample variance, random 
noise and the interdependence of the model parameters. In 
addition to the improvements from this revised analysis, our 
new results also reflect the increased sensitivity of the COBE 
data after two years of observing, along with the more ac- 
curate estimate of the cosmological signal in the Tenerife 
data. 



4.1 Properties of the two data sets 

The instrumental profile of the COBE data is described ap- 
proximately by a Gaussian beam with FWHM~ 7°. As can 
be seen in Figure 1 of Watson et al. (1992) there is a range 
of angular scales to which the COBE DMR and Tenerife ex- 
periments are both sensitive. Measurements taken by COBE 
cover the full sky, but to determine the CMB fluctuations the 
region of the Galactic plane has been excluded (|6| < 20°) 
thereby introducing a degree of uncertainty in estimating 
the properties of the global field; this effect is commonly 
termed sample variance (Scott, Srednicki & White 1994). 
The uncertainties in the COBE two year results are domi- 
nated by the effect of cosmic variance i. e. the fact that our 
stochastic theory describes the Universe as a particular re- 
alisation of a random field. Together the cosmic and sample 
uncertainties form an intrinsic limitation of the COBE ex- 
periment since unlike random errors they are not reduced by 
increased integration time. The two year COBE data have 
been analysed independently by a number of authors (e.g. 
Banday et al. 1994, Bennett et al. 1994, Bond 1994, Gorski 
et al. 1994, Wright et al. 1994). All find evidence for statis- 
tically significant structure at an amplitude consistent with 
that of 30 ± 5 nK rms on a 10° scale announced by Smoot 
et al. (1992) for the first year data. The best fit values for 
n and Qrms-ps depend on the precise analysis techniques 
employed, but are generally consistent with the values of 
n = 1.10±0.29, Qrms-ps = 20.3±4.6 found by Tegmark 
and Bunn (1995) for the combined 53 and 90 GHz data with 
the quadrupole included. In the case of the Tenerife experi- 
ment the double-switching scheme removes the contribution 
of low order multipoles decreasing the cosmic variance of 
the signal on these large scales; the major source of un- 
certainty is produced by the partial sky coverage (sample 
variance) and the instrumental noise. The region observed 
by the Tenerife experiments covers ~ 5000 square degrees 
but here we have limited our analysis to our region of deep- 
est integration at high Galactic latitude which constitutes a 
sample ~ 500 square degrees. For such a region the uncer- 
tainties due to the partial sky coverage dominate over the 
intrinsic variance by a factor ~ 10 (Scott, Srednicki & White 
1994) and the combined uncertainty is approximately of the 
order of that introduced by the instrumental noise in the 
15+33 scan. 

In previous work (Hancock et al. 1994) we explicitly 
took into account the effects of cosmic and sample variance 
by using Monte Carlo simulations. This was necessary be- 
cause the simple excess variance statistic used in the com- 
parison incorporated only the uncertainty due to random 
noise. The comparison by necessity assigned equal probabil- 
ity to all values of n for the COBE data, since the COBE 
first year results were only published in the form of a best 
fit n versus Q rm s- V s relation rather than a full two dimen- 
sional probability distribution. What is required is a data 
analysis technique that allows the joint probability of any 
combination of the model parameters to be calculated and 
which implicitly takes into account random errors and cos- 
mic and sample uncertainties. The Bayesian approach using 
the likelihood function as described in Section 3 attempts 
to do precisely this. The likelihood function peaks at the 
most likely parameters (the best estimate of the true values 
if the likelihood function is unbiased) and has some distri- 
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bution which through Bayes theorem is representative of the 
combined effects of the cosmic, sample and random uncer- 
tainties. The issue of how well this distribution reflects the 
true uncertainties is addressed explicitly in a forthcoming 
paper (Rocha et al. in preparation) by comparison of the 
Bayesian probability distribution with that obtained from 
direct Monte Carlo simulations of the data. The Bayesian 
and frequentist approaches are found to be consistent for 
the Tenerife data and to a good approximation the likeli- 
hood function is also seen to be an unbiased estimator of the 
model parameters. Consequently the likelihood surface for 
the joint Tenerife and COBE data set provides the defini- 
tive means of comparison of the observations under some 
assumed sky model. 

4.2 The Tenerife-COBE likelihood function 

Here we apply the likelihood analysis to the COBE two 
year data and the Tenerife 15+33 scan, assuming a power 
law model with free parameters n and Q rmB - V s- The COBE 
Galaxy-cut two-year map consists of 4038 pixels, whilst the 
Tenerife Galaxy cut (RA 161° — 230°) scan contains 70 pix- 
els, requiring a 4108 x 4108 covariance matrix for a joint 
likelihood analysis of the data. The direct inversion of such a 
large matrix, necessary for the likelihood analysis, is compu- 
tationally intensive, but has been implemented for the com- 
bined 53 and 90 GHz COBE two-year data by Tegmark and 
Bunn (1995), hereafter "TB95". A number of other authors 
(Bond 1994, Bunn and Sugiyama 1994, Gorski 1994) have 
computed the likelihood function by using various data com- 
pression techniques to reduce the size of the covariance ma- 
trix. The compression is achieved by discarding noisy data 
vectors whilst retaining most of the cosmological signal, and 
although the results are close to optimal, the slight loss of 
signal results in marginally larger error bars than the "brute 
force" approach. The latter method is conceptually the sim- 
plest one, since it is merely a likelihood analysis using all 
available data, and involves no adjustable parameters such 
as the degree of data compression. Since each matrix inver- 
sion requires merely about 10 minutes on a fast workstation, 
we use the brute force method here. 

We arrange the pixels in a 4108-dimensional vector and 
compute the likelihood function as in TB95 by Cholesky 
decomposition of the 4108 x 4108 covariance matrix at a 
dense grid of points in the (n, Q r ms- P s)-parameter space, 
marginalizing over the four "nuisance parameters" that de- 
scribe the monopole and dipole. The covariance matrix 
consists of three parts: a 70 x 70 block with the covari- 
ance between the Tenerife pixels, a 4038 x 4038 block with 
the covariance between the COBE pixels, and off-diagonal 
4038 x 70 blocks containing the covariance between the 
Tenerife and COBE pixels. Simply multiplying the likeli- 
hood curves resulting from two separate analyses of the 
COBE and Tenerife data sets would correspond to neglect- 
ing the off-diagonal blocks, and this is clearly only a good 
approximation if the two data sets are almost uncorrelated. 
We find that the inclusion of the cross-terms makes a non- 
negligible difference, which is not surprising in view of the 
fact that the two experiments probe comparable angular 
scales and have observed a common region of the sky. 

For our Tenerife data, the best estimate of the cosmo- 
logical signal is obtained from the 15+33 combined scan 
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Figure 5. The cross-correlation function C(0) for the 15 GHz 
and 33 GHz data over RA=161°— 230° showing the data points 
and the one sigma error-bars. The solid line corresponds to the 
prediction of a Harrison-Zel'dovich spectrum of fluctuations with 
a normalization of Q r ms-ps = 22 fiK. The dashed lines are the 
la confidence bounds arising from cosmic variance and sample 
variance. 

after correction of the error bars for the correlated atmo- 
spheric noise term. The possible contribution of Galactic 
signals has been estimated from the 10 GHz data to be less 
than 4 fiK at 33 GHz and has not been considered in the 
current comparison. The normalised likelihood function for 
this scan, as plotted in Figure 6, represents the joint prob- 
ability of obtaining a given combination of n and Q r ms- P s- 
On its own, the Tenerife configuration provides less lever- 
age on the slope of the spectrum than the COBE satel- 
lite. This is because the Tenerife experiment is insensitive to 
the largest angular scales, and because the one-dimensional 
shape of the dec+40° strip makes it difficult to separate the 
power contributions from different scales. In other words, 
a narrow strip corresponds to wide window functions in t- 
space, with considerable aliasing of small-scale power onto 
larger scales (just as the case is with one-dimensional "pencil 
beam" galaxy surveys). As a result, the Tenerife data can be 
equally well fit by a range of n and Qrms-ps values, result- 
ing in a likelihood ridge in (n, Q rms _ ps )-space with minimal 
discriminatory power for the parameter n. In contrast, the 
COBE observations are sensitive to the slope to the extent 
that the likelihood surface is peaked in the n-dimension. 
Combining the COBE information with the Tenerife data 
improves the situation in two ways: it extends the lever arm 
on the spectral slope from the COBE scales down to the 
4° scale of Tenerife, and in addition eliminates the above- 
mentioned aliasing problem, since the joint data set is no 
longer one-dimensional. 

Figure 6 shows the confidence contours obtained from 
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Figure 6. Constraints on the quadrupole Q r ms—ps and on the 
spectral index n of fluctuations obtained from a joint likelihood 
analysis of the Tenerife data and COBE DMR two-year data. 
The contour levels represent 68 %, 95 % and 99 % of the region 
of joint probability. The peak of the distribution lies at n = 1.37, 
Qrms-ps = 16.1/iK and is identified by the cross. 

Bayesian integration under the combined COBE and Tener- 
ife likelihood surface assuming a uniform prior. The 68% 
joint confidence region in (n, Q rms _ ps )-space encloses a re- 
gion from 0.90 to 1.73 in n for Qrms-ps in the range from 
12.1 to 22.9 fj,K, with the peak at n — 1.37, Qrms-ps = 
16.1/iK. Marginalizing over Qrms-ps with a uniform prior, 
one obtains the probability distribution for n as given in Fig- 
ure 7, corresponding to n = 1.33 ± 0.30 at 68% confidence. 
The resulting limits on the normalization, conditioned on 
n = 1 as is customary, are Q rm s-ps = 21.0 ± 1.6. The cor- 
responding results in TB95 using just the COBE data, and 
including the weak correlated noise term (Lineweaver et al. 
1995), were n = 1.10 ±0.29 and Q r ms- P s = 20.3 ± 1.5, with 
the peak likelihood located at n = 1.15, Qrms-ps = 18.2/^K. 
In other words, although the total normalization has risen 
by a mere 3%, the slope estimate has risen by 12% and 
the peak likelihood has been shifted to higher n and lower 
Qrms-ps- This indicates that the higher angular resolution 
data from the Tenerife experiment contains slightly more 
power on small scales. As explained in Section 3.1, this is 
not unexpected, since the presence of a Doppler peak would 
cause a rise in the power spectrum at higher I. 



5 ADDITIONAL CONSTRAINTS ON THE 
PRIMORDIAL SPECTRAL INDEX 

The large angular scale CMB observations from Tenerife 
and COBE probe fluctuations that have yet to go non-linear 
and the shape of the power spectrum is thus insensitive to 
the exact abundance of the baryonic mass (fi&) and to the 
value of h = Ho/100 km s _1 Mpc~ . However, together the 
Tenerife and COBE observations provide a direct measure 
of the CMB power spectrum normalisation, against which 
one can compare intermediate scale observations and hence 
discriminate between competing cosmological models. Here 
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Figure 7. The marginal likelihood for the spectral index n as 
obtained from the joint analysis of the Tenerife and COBE data. 
The spectral index is seen to lie in the range 1.02 < n < 1.62 at 
68% confidence, with a best fit value of n = 1.33. 

we stress the importance in the agreement between the de- 
rived normalisation for the independent Tenerife and COBE 
experiments, which are subject to different systematic er- 
rors and different foreground contamination. The accuracy 
to which we know the normalisation of the power spectrum 
clearly becomes an issue when comparing with smaller scale 
observations to determine cosmological parameters. This is 
a particular concern if, as has been suggested (Crittenden 
et al. 1993, Steinhardt 1993, Abbott and Wise 1984), a com- 
ponent of the large scale anisotropy may be due to tensor 
metric perturbations produced from a background of gravi- 
tational waves. These can arise naturally in inflationary sce- 
narios and would contribute a component C i ' T ' =< \af m \ 2 > 
to the observed CMB angular power spectrum (cf Equation 

(T) 

. The ratio of the tensor modes C, to the scalar modes 
Cy^' (primordial density fluctuations) is highly suppressed 
for fluctuations contained within the horizon volume at re- 
combination and hence their contribution is only significant 
on scales ~ 2°. Consequently the existence of a tensor con- 
tribution has implications when comparing the large scale 
anisotropy level with that on smaller scales and with large 
scale structure observations in order to test cosmological 
models. The anisotropy measurements from Tenerife and 
COBE fix the sum of C\ and C\ , but the separation of 
the two terms requires a comparison with smaller scale ob- 
servations under the assumption that a given cosmological 
model is correct (Crittenden et al. 1993, Steinhardt 1993). 
In the case of power law inflation a relation exists between 
the tensor to scalar ratio /Cj an d the slope of the pri- 
mordial power spectrum (Crittenden et al. 1993, Steinhardt 
1993) : 

Cf'/Cf' «7(l-n) (4) 
from which we see that the two contributions are comparable 

(T) (S) 

at n = 0.85 with C\ decreasing relative to Cy for higher 
values of n. Thus the limit of n > 1.0 obtained from the 
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ana lysis of the combined Tenerife and COBE data in Section 



4.2 implies that it is unlikely that the tensor component will 
be dominant in such models. In the remainder of this section 
we shall investigate this in more detail by comparing with 
medium-scale anisotropy results. 

Contemporary cosmological models with adiabatic fluc- 
tuations predict a sequence of peaks on the power spectrum 
which are generated by acoustic oscillations of the photon- 
baryon fluid at recombination. Of particular interest is the 
height and position of the main acoustic peak — the so 
called Doppler peak: the height depends on quantities like 
the baryonic content of the universe, fi&, and Hubble con- 
stant, Ho, whilst the position depends on the total density 
of the Universe, Slo . 

The work required to place the existing medium scale 
anisotropy results into a common statistical framework and 
then to compare them with the predictions of cosmological 
models, has been carried out by Hancock et al. (submitted) 
who find strong evidence for the existence of a Doppler peak 
on medium scales. Further details together with a compar- 
ison with a fuller range of cosmological models will be pre- 
sented in Rocha et al (in preparation). For the current pur- 
pose of discussion of the value of n, and comparison with the 
results obtained from Tenerife and COBE alone, we present 
here a version of the results for n found in Hancock et al. 
(submitted) . 

The precise form of the Doppler peak depends on the 
nature of the dark matter, and the values of Qo, and 
Ho. Thus in order to use the medium-scale anisotropy re- 
sults to provide additional leverage on the spectral index 
determination it is necessary to adopt a given cosmological 
model. Taking the minimalist assumption of the standard 
Cold Dark Matter (CDM) model, without gravity waves, 
with Ho = 50kms _1 Mpc -1 , Q,b = 0.07, Hancock et al. (sub- 
mitted) find this offers a good fit to all of the data. Tilting 
the initial spectrum of fluctuations away from scale invari- 
ance and fitting to the ATi we can delimit n using the full 
data set. 

Hancock et al. considered the four-year COBE data, and 
used the binned angular power spectrum (Tegmark 1996) in 
conjunction with Tenerife, Python, South-Pole, Saskatoon, 
MAX, ARGO, MS AM and CAT experiments. As seen in 
Figure 8, only models with initial spectra in the range 1.0 < 
n < 1.2 are allowed by this improved analysis, with best 
fit values of Q rm s-ps = 15 [iK for n = 1.1. These results 
are consistent with those from large scales alone, and with 
the inflationary value of n ~ 1.0. We note that the presence 
of a gravity wave component in our model would require 
even larger values of n than those derived above and given 
Equation ^ this reaffirms our conclusion that a significant 
gravity wave contribution is unlikely. The above constraint 
on n depends on which values of fio, ^6 and Ho are chosen, 
but is not simply related to any one parameter. Constraints 
on n derived using a fuller set of models are given in Rocha 
et al (in preparation). 

Prior to the COBE detection the most common nor- 
malization of models used was based on the value of a%{ 
(AM/M) r ms in a sphere of radius 8ft _1 Mpc) derived from 
galaxy clustering assuming a bias b — 1 i.e. that light traces 
mass (Kaiser 1984, Davis et al. 1985). If in fact galaxies are 
more highly clustered than matter, the amplitude of the ini- 
tial matter perturbations (and hence AT rms /T ) necessary 
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Figure 8. Recent CMB anisotropy results on large and interme- 
diate angular scales are used to delimit the spectral index n of 
the primordial fluctuations. 



to produce the observed clustering is reduced by the factor 
b. The Tenerife and COBE data provide an independent and 
more accurate value for the normalisation which for a given 
cosmological model allows us to determine the degree of bias 
necessary for consistency with the large scale structure ob- 
servations. 

In particular the COBE normalized CDM model with 
a tilt of n = 1.1 gives a bias of the order ~ 0.7. 



6 DISCUSSION AND CONCLUSIONS 

We have presented a detailed analysis of the measurements 
taken at Dec=+40° by the Tenerife experiments at 10, 15 
and 33 GHz. After accounting for both local atmospheric 
and discrete radio source foregrounds, the 15 and 33 GHz 
data at high Galactic latitude are seen to contain statis- 
tically significant signals which have their origin in com- 
mon hot and cold features. The cross-correlation function 
between the data at these two frequencies demonstrates 
that the amplitudes and shapes of the structures detected 
at 15 and 33 GHz are similar. This, combined with our 
measurements at the lower frequency of 10 GHz, implies 
that the CMB signal dominates over the Galactic contribu- 
tion at 15 GHz, and that the maximum possible Galactic 
contribution at 33 GHz is smaller than 10 % of the de- 
tected signal. Our best estimate of the cosmological signal 



is Qrms- P s = 22_ 6 fiK for an n = 1 inflationary spectrum. 



This amplitude is reduced by 4fiK over that previously re- 
ported for the same data set and results from an improved 
separation of signal from atmospheric noise. Comparison 
of the Tenerife and COBE two year anisotropy detections 
by means of the likelihood function allows a detailed in- 
vestigation of the allowed parameter space for a power law 
model of the fluctuation spectrum. The best fit values of 
n and Qrms-ps are 1.37 and 16/iK, and marginalising over 
Qrms-ps we find both data sets consistent with n in the 
range 1.0 < n < 1.6. These results support inflationary 
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models, which predict n ~ 1 and future improvements in 
the Tenerife data and use of the 4 year COBE data will 
narrow the range of allowed n. Improvements of this kind 
are important to determine the power spectrum normalisa- 
tion, since only on these large angular scales is it possible 
to place limits on the tensor to scalar ratio independent of 
the precise details of the cosmological model. By combining 
the large scale anisotropy measurements from Tenerife and 
COBE with observations of medium scale anisotropy mea- 
surements, improved limits have been placed on n, but at 
the expense of assuming an underlying cosmological model 
(in our case CDM). 

For a CDM model with H = SOkms" 1 Mpc~\ fl b = 
0.07, we find that this data set is consistent with n in the 
range 1.0 < n < 1.2. The best fit values of n and Q rm s- V s are 
1.1 and 15/iK. This COBE normalized tilted model predicts 
a bias of the order ~ 0.7. So although this model fits well the 
CMB data alone it may not give a realistic scenario when 
we consider jointly the CMB data and the observed galaxy 
clustering. 

Our attempts to place increasingly more accurate lim- 
its on fundamental cosmological parameters will undoubt- 
edly place increasing demands on observers for ever im- 
proved accuracy until the point where the intrinsic cos- 
mic variance becomes the dominant form of error. Confi- 
dence in the results to this tolerance level will probably re- 
quire close co-operation between observers and the combina- 
tion of results from space-based, balloon-based and ground- 
based telescopes working over a range of frequencies. Poten- 
tially the most powerful observations will result from map- 
ping overlapping/ interlocking CMB fields with independent 
multi-frequency instruments. In contrast to the statistical 
results currently being reported, this method builds in the 
necessary redundancy to reduce systematic errors to an ac- 
ceptable level. Such observations are currently in progress 
with the Tenerife instruments, which in conjunction with the 
COBE 4 year data should provide a high signal to noise map 
of the last scattering surface thus providing a two dimen- 
sional representation of the seed structures as compared to 
the simple 1-D scans reported here. Being on scales greater 
than the horizon size at recombination the form of such a 
map would reflect the structures generated in an inflationary 
driven phase in the very early universe. 
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